Lowest Landau level broadened by a Gaussian random potential with an arbitrary 
correlation length: An efficient continued-fraction approach 
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For an electron in the plane subjected to a perpendicular constant magnetic field and a homoge- 
neous Gaussian random potential with a Gaussian covariance function we approximate the averaged 
density of states restricted to the lowest Landau level. To this end, we extrapolate the first 9 coeffi- 
cients of the underlying continued fraction consistently with the coefficients' high-order asymptotics. 

• We thus achieve the first reliable extension of Wegner's exact result [Z. Phys. B 51, 279 (1983)] for 
' the delta-correlated case to the physically more relevant case of a non-zero correlation length. 
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Nearly ideal two-dimensional electronic structures 
' have attracted great attention for more than a decade 
, not only because of their varied and important applica- 
' tions, but also because of the discovery of the quantum 
' Hall effect Q. For a microscopic understanding of the 
, ' occurring phenomena it is essential to know the spectral 
, properties of electrons confined to two dimensions under 
the influence of a perpendicular constant magnetic field 
^ taking into account the presence of disorder. 

A commonly studied minimal model is that of non- 
interacting electrons which is characterized by the one- 

• electron Hamiltonian given by the Schrodinger operator 

0\ '. 

H:=K + V (1) 
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Here x :— (xi^X2) are Cartesian coordinates of the Eu- 
clidean plane M^, h is Planck's constant, e is the elemen- 
tary charge, m is the (effective) mass of the (spinless) 
electron, and B > the strength of a perpendicular mag- 
netic field. The static random potential V is added to 
mimic the interaction with quenched disorder. Its prob- 
ability distribution is assumed to be Gaussian with zero 
mean and Gaussian covariance, that is 



V{x) = 0, V{x)V{x') = cr^ exp 



2A2 



(3) 



The overbar denotes averaging with respect to the prob- 
ability distribution. Finally, ct > is the strength and 
A > the correlation length of the fluctuations of the 
potential. 

The spectral resolution of the unperturbed Hamilto- 
nian K reads K = X^J^o^"-^"' "^tiere the eigenvalue 
e„ := (2n + 1)-^^ is called the n-th Landau level and En 
denotes the corresponding eigenprojector. With increas- 
ing B the magnetic length I :— ^Jh/{eB) decreases, while 
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the degeneracy (a;|£'„|a;) = {2itI'^)~^ (per area) of each 
Landau level and the distance between successive levels 
increase. Hence, for a fixed concentration of electrons, 
sufficiently high fields, and low temperatures it is reason- 
able to simplify the model by restricting its Hamiltonian 
(|I|) to the (still infinite dimensional) eigenspace EoL'^(R'^) 
of the lowest Landau level. This has been done, for ex- 
ample, in the important work of Wegner Q, where the 
averaged restricted density of states (per area) 



g{e) {y\EoS{e - E^HE^)Ea\y) (4) 

has been calculated exactly in the delta-correlated limit 
A I 0, (T ^ cx), A^CT^ = const. 

The purpose of the present Letter is twofold. One goal 
is to present what we think is an accurate extension of 
Wegner's result to arbitrary values of a and A. Clearly, 
this is of physical interest because in the high-field limit 
the actual correlation length is no longer small in com- 
parison with the magnetic length. The other goal is to 
demonstrate thereby that the power of the well-known 
continued-fraction approach to spectral densities of non- 
trivial quantum problems can considerably be en- 
hanced in cases, where one has an a-priori information 
about their asymptotic high-frequency decay. 

To get rid of physical dimensions we write g in the 
standardized form 



(5) 



In this way we have introduced a one-parameter family 
of even probability densities on the real line R with nor- 
malized second moment 

W{a,u) = W(a,-u) >0 (6) 



duW{a,u)^l= duW{a,u)u^. (7) 



Here we have used Eq. (7) of Q to normalize the second 
moment. 

In this notation Wegner's result for the delta-correlated 
limit reads [0161 



l+[27r-V2/;deexp(e2)]' 

For the other extreme of the spatial extent of correlations, 
namely the constantly correlated case A = oo, one simply 
has the Gaussian 

W{oo, u) = (27r)"i/2 exp(-uV2) • (9) 

For intermediate values of a no exact expression for 
W{a,u) is known. What is exactly known, however, is 
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the fact that W(a, u) falls off for sufficiently large |u| like 
a Gaussian. More precisely, by Eqs. (17) and (15) of 
(see also |^) 

lim ^\n{W{a,u)) = -^±^. (10) 
ti— >±oo tt^ 2a + 2 

In the sequel we will design approximations to the 
Stieltjes transform 



R(a, z) := [ du — ^ \ \ Rez > 0, 
Jm z-iu 



(11) 



of the standardized density of states W which in turn 
will yield approximations to W by means of the inversion 
formula 

W{a,u)^-Vm\Re[R{a,v + m)]. (12) 

According to Stieltjes' classical theory, see for example 
|^,|o|, R can be expanded into a Jacobi-type continued 
fraction 

i?(a,z)= K(^), r,(a)>0. (13) 

J — z 

Here we are using the notation 
oo A . 1 

A'T> - A. <») 
z H A 



z+ ■■ 

for the continued fraction with coefhcients Ai,A2,... 
and variable z. 

To derive the continued-fraction coefficients {rj{Q)} 
and {rj(oo)} corresponding to (H) and (||), respectively, 
we use the identity 

K(jI±2i), ^;f/^'--f:;;;> :r(ft,.,). as) 

vahd if 7 > 0, /3 + 7 > 0, and Rez > 0. Here Z?^ denotes 
Whittaker's parabolic cylinder function with index see 
Section 8.1 of jll]. The identity follows by iteration from 
the observation that the rhs of (ISf) obeys 



r(/3, 7, z) = [z + (/? + 7)r(/3 + 7, 7, z)]-i (16) 

which itself is a consequence of the first of the four re- 
currence relations for the Z),y's in Section 8.1.3 of ||ll[. In 
fact, Eq. ( p^ is equivalent to Eq. (14) in §50 of [§. With 
the help ofQ one now checks that 

^j(0)-^ + ^i, r,(oo)=j, (17) 
for all j > 1 . 
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According to ||T^ the asymptotic behavior ( |lO| ) imphes 
the following asymptotic linear growth for the continued- 
fraction coefhcients 

To (a) a + 1 , ^ 

lim = ——. (18) 

j^oo J a + 2 

It is natural to match this linear growth to the first 
J < oo coefficients ri{a), . . . , rj(a) to construct the an- 
nounced approximations R^''^\a,z) to R{a,z) by means 
of 

R(-')^a,z):^ K {^^), (19) 
j=i z 



where 



(J), X / rj{a) for j < J 



'^^^'=[rAa) + mU~J) for ^''^ 

Not surprisingly, the approximation i?*^'^' to R results in 
an approximation W'^-'^ to W with the same Gaussian 
fall-off ( |lO|) in the tails. This can be seen as follows. By 
virtue of (ij) the approximation can be expressed 
as a terminating continued fraction 

RiJ)(a,z) = 



z+ ''^"^ 



z + rjia)Tirj{a),^,z) 

(21) 

In view of (O) it is therefore sufficient to show 



and 



liin ±ln{Re[T{f3,j,iu)])^~l- (22) 

u—^±00 U 2^ 



lim Im[r(/3, 7, iu)] 0, (23) 

u — >±oo 

because these asymptotic properties are conserved under 
the (repeated) substitution 

T(/3, 7, iu) ^ [iu + aT{(3, 7, iu)]-^ (24) 



for any a > 0. The validity of (|2^) und ( p3| ) can be de- 
duced from an asymptotic evaluation of the following 
Riccati differential equation 

T^T'l/S, 7, iu) = i/3 [T(/3, 7, m)]^ - uT(/3, 7, iu) - i 

(25) 

which itself follows from the aforementioned recurrence 
relations. To summarize, Eq. (15) adds a two-parameter 
family of terminators for affine-linear extrapolations to 
the toolkit of 01. 
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In order to compute the first J continued-fraction coef- 
ficients ri(a), . . . , rj(a), we employ the well-known one- 
to-one correspondence [p|,p^P] to the first J even mo- 
ments M2(a), . . . , M2j{a) of the standardized density of 
states. 

According to (|^) and (||) one has for the latter 
M2j{\^/f) -.^ [ duW{\^/f,u)u'^^ 



2 /\2\3 



= 2nl'i^ {y\{EQVE:,YJ\y) (26) 

if j > 1. We now substitute V = J^2 x V {x)\x) {x\ into 
the rhs of (|2q) and use the standard reduction formula 



2] (2j-l)!! 3 

nmfc))- E Y[^i<Pk{2s-mvix{Pki2.s))) 

fc=l k=l s=l 

(27) 

for the average of a product of 2j jointly Gaus- 
sian random variables with zero mean. Here 
x{l), . . . ,x{2j) G are 2j points of the plane and 
Pk denotes the fc-th of the (2j — 1)!! permutations of 
the first 2j natural numbers which lead to different 
sets {{Pfe(l),Pfe(2)}, . . .,{Pk{2j - l),Pfc(2j)}} oij pairs 
{Pfe(2s — 1), Pfc(2s)}. Since the covariance function ^ of 
the random potential V and the position representation 



^(a;i4 - X2x[) - -^{x - x')' 



(28) 



of the eigenprojector corresponding to the lowest Landau 
level are both Gaussian, one ends up with a sum over 4j- 
dimensional Gaussian integrals which can be performed 
to yield 

/ 1 
M2j(a)= 1 + - > ^ . , , , a>0. (29) 
^ V aJ ^ det(Afe(a))' 

Here the 2j x 2j-matrix A^ia) is defined through its en- 
tries in terms of the Kronecker delta 



1 1 

(^fe(a))p,^ := (l + - ^A'+l.i' - -E ('^A',Pfc(2s-l)VPfc(2s) + Si.^P^(2s-l)S^,PU2s)) 



s=l 



(30) 

We have computed the sum (|2£ ) for j = 1, 2, . . . , 6 us- 



ing the Axiom symbolic system [1^ for general a. The 
results for M2j{a), j = 1, 2,3,4, are given by the follow- 
ing rational functions 
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M2(6- 1) = 1 
M4(5-l) = (362 + 2)/(62 + i) 
Me{b - 1) = {15b^ + 75b* + 1026^ + 30)/ {b^ + 6b* + llb^ + 6) 

Ms{b - 1) = (105624 + 2835622 + 3283062° + 2136976i8 + 8611306i6 + 22318076i4 
+3750338612 _,_ 40385436i° + 27171176® + 11055106^ + 2612166^ + 3279262 + 1680)/ 
(624 _,_ 29522 ^ 355^20 _,_ 262061® + 1185461^ + 352766i4 + 699746^2 + 919066i° 
+780256® + 410156^ + 124616* + 195462 + 120). 

(31) 

Since the expressions for Mio(a) and Mi2{a) are ex- 
tremely lengthy, we only give, as an example, their values 
for a = 1: 



Mio(l) 



158 659 605 940 126 452 841 



294 310 802 651335 470 

n,r /-|^ _ 388 336 271 072 847 928 549 926 597 113 071 401 088 677 997 478 405 727 555 223 031 
-""J 12 1, 1 j g2 363 680 790 265 452 914 044 225 729 941 466 953 642 191 484 142 275 602 750 • 

(32) 

The moments Alicia) and MiQ{a) have been computed 
exactly as a reduced fraction for a — j,^, 1,2,4 only, by 
employing a C-program. For the same set of values for a 
we also have computed Mis{a), but only have been able 
to gather up the 17!! « 3.4-10'' terms occurring in ( p9| ) to 
a sum of approximately 10* reduced fractions, the value 
of which has been computed accurately by floating-point 
arithmetics. For lack of space we only show the first 35 
digits for the example a = 1: 

Mi4(l) = 47657.946072630475536554559639613945 . . . 
Mi6(l) = 545841.43592501744224295351679205436. . . 
Mi8(l) = 6980770.3795705620571551127542699202 . . . 

From Fig. ^ it can be seen that the resulting 
continued- fraction coefficients ri (a) , . . . , rg (a) approach 



quite rapidly their asymptotic behavior given in (18). 
For example, the relative extrapolation error |rg (a) — 
rg{a)\/rr){a) varies between 0.2% for a = 4 and 0.02% for 

In Fig. we demonstrate the convergence of the result- 
ing approximations M^*^-^-* (1, u) to W{l,u) for increasing 
J — 1,2,. ..,9. The differences of two successive approxi- 
mations, W'^'^''{a,u)—W^'^~*'>{a,u), seem to form a nearly 
alternating sequence with geometric decrease for fixed u. 
Both observations, taken together, suggest rapid point- 
wise convergence of the sequence {W'^J)}. Therefore, we 
may safely conclude that T4^(^) constitutes a reliable ap- 
proximation to W. Figure p| shows a plot of the corre- 
sponding approximation g^^> to the averaged density of 
states g for different A. Note that, by construction, p*-^-* 
is exact in the limiting cases of a delta-correlated and 
a constantly correlated random potential; more impor- 
tantly, e ^ 277^2 p(9)(e + eo) is an even probability density 
with the same first 18 moments and the same Gaussian 
fall-off in the tails as the true density for general values 
of the correlation length A. 
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FIG. 1. Symbols represent the first 9 continued-fraction 
coefficients ri(a), . . . , rc){a) for different values of a. Note that 
ri (a) = 1 for all a by (Q) The straight lines correspond to the 
asymptotics of {rf\a)} as defined in (EOh. 
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FIG. 2. Difference of two successive approximations to 
for a = 1 as a function of w given on a linear (a) and a 
logarithmic (b) scale. 
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FIG. 3. The approximation ^'^^ to the averaged density 
of states Q for different values of the correlation length A as a 
function of the energy e. 
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